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A review of spectroscopic results obtained from Chandra High Energy Transmission Grating Spectrometer 
(HETGS) and XMM-Newton Reflection Grating Spectrometers (RGS) observations of several wind-fed high-mass 
X-ray binaries (HMXBs) is presented. These observations allow us to study the structure of the stellar wind in 
more detail and provide, for the first time, a dyanmical view of the X-ray photoionized wind that surrounds the 
compact object. At the same time, however, they are also providing us with numerous puzzles that cannot be 
explained in terms of simple models. For example, simple spherically-symmetric wind models cannot explain the 
observed orbital-phase variability of the line intensities and shapes, which may be caused by intrinsic asymmetries 
due to the presence of the compact object and/or more complicated radiative transfer effects. The observed line 
shifts are smaller than those expected from extensions of simple wind models of isolated OB supergiants. In 
addition, several novel spectroscopic discoveries have been made, including: (1) P-Cygni lines from an expanding 
wind, (2) detection of multiple Si K fluorescent lines from a wide range of charge states, (3) Compton scattered 
Fe K lines from a cold medium. We discuss how these spectroscopic diagnostics can be used to understand some 
of the global properties of stellar winds in HMXBs. 



1. Introduction 

In wind-driven high mass X-ray binaries 
(HMXBs), a neutron star or a black hole sweeps 
up a small fraction of a stellar wind lost by a 
massive O- or B-type companion star. As mat- 
ter is accreted onto the compact object, a frac- 
tion of the gravitational potential energy is con- 
verted into X-rays, which then ionizes and heats 
the surrounding gas. The wind reprocesses hard 
X-rays from the compact object, resulting in dis- 
crete emission lines and continuum radiation that 
carries a wealth of information about the physical 
state of the reprocessing medium. The compact 
object can, therefore, be used as an illuminating 
source to probe the structure of the stellar wind 
and derive the physical parameters that charac- 
terize its nature. 

In addition to heating and ionizing the stel- 
lar wind, the X-ray source is also responsible for 
dynamically affecting the wind through its in- 
tense gravitational field and by destroying ions 
with strong UV resonance transitions that drive 
the stellar outflow. These effects, in combination 
with the orbital motion of the binary system, also 
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produce a "wake" that trails behind the compact 
object. In low-luminosity HMXBs, however, the 
degree of stellar wind disruption by the compact 
object is expected to produce only a minor per- 
turbation to the observed X-ray spectrum. This 
is especially true when observed during eclipse, 
since many of the dynamical effects are produced 
primarily in the vicinity of the compact object. 

ASCA observations of several HMXBs have 
shown that their X-ray spectra exhibit both soft 
X-ray emission from highly ionized ions and flu- 
orescent lines from cold, less ionized material 
(Vela X-1 - Nagase et al. 1994; Cen X-3 - Ebi- 
sawa et al. 1996; GX301-2 - Saraswat et al. 
1996). Although the dominant excitation mech- 
anism (i.e., coUisional or photoionization-driven) 
that is responsible for producing the soft X-ray 
lines cannot be determined unambiguously from 
these observations, cascades following recombi- 
nation seemed to be the most natural candidate 
given the presence of an intense X-ray continuum 
radiation field. Subsequently, Liedahl & Paerels 
(1996) and Kawashima & Kitamoto (1996), for 
the first time, detected a narrow radiative recom- 
bination continuum (RRC) of S XVI in the ASCA 
spectrum of Cyg X-3, which provided unequivocal 
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Figure 1. A schematic drawing showing the 
various componets of an intermediate-luminosity 
{Lx ^ 10^^ erg s'^) HMXB system. For more lu- 



minous systems, the dense clumps may also be 
highly ionized, emitting H- and He-like recombi- 
nation lines and showing very little fluorescence 
emission. 



evidence that the lines from highly-ionized gas in 
Cyg X-3, and most likely in all HMXBs, are driven 
almost exclusively through photoionization. 

The first physically-motivated modeling of the 
X-ray spectrum was presented by Sako et al. 
(1999) using the same ASCA spectrum of Vela X- 
1 originally published by Nagase et al. (1994). 
Sako et al. (1999) have argued that the presence 
of emission lines from both high-ionized and cold 
near-neutral material can be understood only if 
the wind is structurally inhomogeneous, consist- 
ing of cool dense clouds embedded in a hot, highly 
ionized medium. They characterized the wind 
velocity profile as derived by Castor, Abbott, & 
Klein (1975; hereafter, CAK), 
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where i?* is the radius of the companion star and 
Voo is the terminal velocity and /3 is a number that 
typically Hes in the range 0.5 — 1.0 (Abbott 1986; 
Pauldrach, Puis, & Kudritzki 1986). The particle 
density is, then, uniquely specified everywhere in 



n{r) 



M 



1 



r 



(2) 



the wind for an assumed 
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and, hence, the ionization parameter ^ = 
Lx/inr^) is also uniquely determined for a given 
X-ray luminosity. Here, r is measured from the 
center of the companion star and /i is the aver- 
age atomic mass unit of the constituent particles 
(/J, ~ 1.3 for solar abundances). Using this rela- 
tively simple model for the wind dynamics, they 
were able to infer the velocity structure of the hot 
medium and the mass-loss rate associated with 
that component. They also estimated the total 
mass in the clouds through measurements of flu- 
orescent line intensities, and concluded that the 
cold fluorescing material contains a large fraction 
of the total wind mass (> 90%), while most of 
the wind volume (> 95%) is occupied by highly 
ionized matter. According to this picture, the 
stellar wind in Vela X-1 is driven by UV radia- 
tion from the star that provides a driving force 
on the clumpy material, while the hot, ionized 
component is essentially transparent to UV radi- 
ation. Subsequently, Wojdowski, Liedahl, & Sako 
(2001) presented a detailed analysis of the X-ray 
spectrum of Cen X-3, which is on average more 
luminous than Vela X-1 by roughly an order of 
magnitude, and showed that the dense clumps 
that produce fluorescence lines in Vela X-l-like 
systems are more highly ionized in Cen X-3. The 
mass loss rate inferred from H- and He-like ions 
alone was consistent with those of normal isolated 
O stars. 

While the most general characteristics of stel- 
lar winds in HMXBs were revealed using moder- 
ate spectral resolution data acquired with ASCA.^ 
there are still many outstanding issues that can 
only be addressed with high resolution spectro- 
scopic data available from Chandra and XMM- 
Newton. The most significant improvement per- 
haps is the ability to study the dynamics of 
the X-ray emitting gas, which cannot be stud- 
ied with the ASCA detectors, because line shifts 
and widths due to typical wind velocities of ~ 
1000 km s~^ require a spectral resolving power 
an order of magnitude higher than that of CCDs 
(i? of at least ~ 300). Accurate column density 
measurements are also possible through measure- 
ments of line ratios of the He-like triplets and 
other transitions (Kinkhabwala et al. 2002). In 
the following sections, we discuss some of the re- 
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cent advances made with observations of HMXBs 
with the grating spectrometers on Chandra and 
XMM-Newton. 

2. Some General Spectral Properties of In- 
dividual Sources 

Shown in Figure 2 are the Chandra HETGS 
spectra of four objects Cen X-3, Cyg X-3, Vela X- 
1, and GX301-2. As immediately evident from the 
figure, the general properties of both the contin- 
uum and line emission are vastly different among 
each of the sources. The spectrum of Cen X-3 
(1st panel in Figure 2), which is a relatively high- 
lumonisity system {Lx ~ 10"^* erg s~^) contain- 
ing a pulsar with a O 6-8 III companion, is com- 
pletely dominated by continuum emission with a 
relatively low line-of-sight column density of cold 
material (A^h ^ 10^^ cm~^). The spectrum ex- 
hibits very weak absorption lines mostly from H- 
like ions with column densities of ^ 10^^ cm~^. 
During eclipse, however, the spectrum is dom- 
inated almost entirely by line emission as pre- 
sented by Wojdowski et al. (2003). 

In Cyg X-3 (second panel in Figure 2), the 
emission line equivalent widths are larger (see, 
Paerels et al. 2000). This suggests that a sub- 
stantial fraction of the continuum radiation is re- 
processed by the wind, which implies that the 
product of the covering fraction and the column 
density of highly ionized material (Af2 x TVj) is 
higher. Since the X-ray luminosities of Cyg X-3 
and of Cen X-3 during these observations were 
similar, the observed differences in the line spec- 
trum also suggest that the average density of the 
X-ray emission line regions in Cyg X-3 is higher. 
This is in qualitative agreement with the high 
mass-loss rate of the Wolf-Rayet companion and 
the compacteness of the binary system (van Kerk- 
wijk 1995). 

Local anisotropics in the distribution of absorb- 
ing material around the X-ray source can also af- 
fect the line equivalent widths. In the Vela X-1 
spectrum (third panel in Figure 2), the observed 
continuum is highly absorbed while the lines long- 
ward of ~ 6A have high equivalent widths. This 
implies that the emission line regions arc illumi- 
nated by a continuum that is less absorbed than 
what we observe. In other words, the absorber 
that covers the continuum source along our line 
of sight does not block the emission line regions. 
As in Cen X-3, the spectrum of Vela X-1 during 



eclipse is also dominated by emission lines pro- 
duced through photoionization of the extended 
stellar wind, as shown by Schulz et al. (2002). 

The soft X-ray photons of GX301-2 (last panel 
in Figure 2) suffers extremely high intrinsic at- 
tcntuation (Nh ^ 10^^ cm^^). The weaknesses 
of the soft X-ray lines here suggest that both 
the line and continuum emission regions are ab- 
sorbed, most likely by the same material. Also 
note the differences in the iron Ka line equivalents 
widths. The more luminous systems (Cen X- 
3and Cyg X-3) with Lx ^ 10^*^ erg s^^ show 
much weaker lines compared to the less luminous 
systems (Vela X-1 and GX301-2), which have 
Lx ~ 10^7 erg s^^. 
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Figure 2. Global spectral properties of four 

HMXBs observed with the Chandra HETGS. Note 
the differences in the continuum and emission line 
properties. 



The X-ray light curves differ substantially 
among the various systems as well. Many sources 
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Figure 3. X-ray light curves of Vela X-1, an X-ray 
pulsar with a period of ~ 280 sec (top), 4U1700- 
377, a neutron star system that shows no pulsa- 
tions (middle), and Cyg X-3, a black hole candi- 
date (bottom). 



Figure 4. Observed Si XIV line profiles in four 
different sources with varying relative contribu- 
tions from emission and absorption. The vertical 
dashed line indicates the rest wavelength of the 
Si XIV line. 



contain pulsars, which show periodic modulations 
of with pulsed fractions of ^ 50% and higher. 
These systems also show flare-like events during 
which the luminosity variabilies by more than a 
factor of few on timescales as short as tens of sec- 
onds, mostly due to accretion instabilites as well 
as variable absorption on longer timescales (see 
Figure 3). On the other hand, there are systems 
like Cyg X-3, a black hole candidate, which gen- 
erally exhibits a smooth X-ray light curve. The 
accretion behavior, therefore, appear to be signif- 
icantly different from source to source. 

3. Wind Dynamics 

One of the most important improvements made 
by observations of HMXBs with Chandra and 
XMM-Newton is the ability to measure velocities 



down to a fraction of the terminal wind velocity. 
In many cases, the observed lines are resolved and 
show complex detailed structure as well as tem- 
poral behavior. 

3.1. Line Profiles and P-Cygni Lines 

The observed line profiles come in many differ- 
ent types as shown in Figure 4, and they provide 
valuable information not only about the dynam- 
ical properties of the emission line regions, but 
they are also sometimes useful for determining the 
dominant excitation mechanism as well. The top 
panel shows a pure emission profile as observed in 
Vela X-1. On the panel below is a P-Cygni line, 
which shows a broad emission feature with a weak 
absorption trough towards the blue. The emission 
component is much stronger than the absorption 
component, suggesting that recombination is the 
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Figure 5. The Chandra data of Vela X-1 observed 
during three different orbital phases (centered at 
(j) 0.0, 0.25, and 0.5; top) and the XMM-Newton 
spectrum observed at = 0.75 (bottom). 



dominant line formation mechanism, since there 
is not enough continuum radiation to account for 
all the observed line photons via resonant fluo- 
rescence scattering. In the third paneP, however, 
the emission and absorption equivalent widths 
are comparable, which suggests that the emission 
component probably comes mostly from scatter- 
ing of the continuum. This interpretation, how- 
ever, can be ambiguous since similar profiles can 
also be generated in a recombination-dominated 
medium with a small covering fraction, thereby 
reducing only the emission component. The last 
panel, finally, shows a pure absorption line profile 
with very little re-emission. In this case, the ab- 
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^Cir X-1 is classified as a LMXB. We, however, show the 
spectrum here since it is the only LMXB system known 
that shows strong P-Cygni lines. 
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Figure 6. Comparison of the global spectra of 
Vela X-1 at orbital phases 4> — 0.0 (red) and 
(j) = 0.5 (black). The fluxes are rescaled to match 
the peaks of the Si XIII forbidden lines of the two 
phases. Note the striking similarity in the emis- 
sion line intensity ratios between the two phases. 



sorber, which happens to lie along the observers' 
line-of-sight, subtends only a small fraction of the 
sky as viewed from the X-ray continuum source. 

3.2. Orbital-phase Variability 

Vela X-1 was observed during three different 
orbital phases centered on = 0.0, 0.25, and 0.50 
with the Chandra HETGS and during (f) = 0.75 
with XMM-Newton. The HETGS spectra shown 
in Figure 5 are easily compared with one another. 
The RGS spectrum is shown separately for clarity. 

A comparison of the observed line spectra at 
(j> = 0.0 and 0.5 is particularly interesting. Shown 
in Figure 6 are the spectra normalized to the peak 
of the Si XIII forbidden line and overlaid on top 
of each other. Apart from the velocity shifts (see 
Figure 7), the spectra look nearly identical with 
^20% difference in the fluxes of most of the lines. 
If anything, the 4> = 0.5 spectrum appears more 
highly-ionized as the Si XIV Lya line is stronger 
and the Si near-neutral line complex is weaker. 
This is in direct contradiction to the spherically- 
symmetric wind model presented in Sako et al. 
(1999), who predicted that the differential emis- 
sion measure distribution at (f> — 0.5 would be 
much softer (i.e., lower ionization) than during 
eclipse due to photoionization of more dense ma- 



6 

terial near the stellar photosphere. The Chandra 
data instead show that there is a large difference 
only in the overall normalization of the distribu- 
tion by roughly an order of magnitude with very 
little difference in the actual shape. However, this 
is not a surprising result, as disruption of the stel- 
lar wind is expected to be most significant near 
the compact object, where both the gravitational 
and X-ray radiation field are high. A simple iso- 
lated wind model, therefore, are not applicable in 
these regions. 
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Figure 7. Line profiles of Lya transitions ob- 
served in Vela X-1 at two orbital phases; 4> = 0.0 
(red) and </> — 0.5 (blue). The fluxes of each line 
pairs are rescaled to match at their peaks. The 
width of the lines are due mostly to the instru- 
ment response. 



Geometrically, the fact that the emission line 
ratios at = 0.0 and 0.5 are similar implies that 
the spectrum emitted in the region blocked by 
the companion during eclipse (a cylinder roughly 
along the line of centers) is nearly identical to 
that of the rest of the wind. A smooth wind 
model characterized by a CAK velocity profile is 
ruled out for reasons stated above. Several simple 
models for the clump distribution can be ruled 
out immediately as well. First of all, a popu- 
lation of identical clumps distributed uniformly 
throughout the wind cannot explain this behav- 
ior, since clumps near the compact object will 
be more highly ionized. In addition, a distribu- 
tion of clumps, whose densities are proportional 
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Figure 8. Si XIV line profiles at observed at 
four different orbital phase ranges in Cyg X- 
3. A particularly interesting feature shown here 
is the difference in the profiles between phases 
(j) — 0.75, which shows a broad, smooth emis- 
sion peak, and = 0.25 that shows what appears 
to be a centrally-reversed line. These differences 
might be due to radiative transfer effects and/or 
geometric occultations by, for example, the com- 
panion star. In either case, the wind must be 
asymmetric with respect to the line of centers. 



to the surrounding local wind density (i.e., higher 
density clumps closer to the stellar surface; see 
Equ. 2), cannot reproduce the data either. In this 
case, we would expect to detect higher line fiuxes 
from lower ionization species during = 0.5, con- 
trary to what is observed. 

Assuming that the wind in the shadow cone 
does not emit a significant amount of X-ray line 
radiation, one expects that most of the wind as 
seen at ^ = 0.0 would be moving away from the 
observer emitting lines shifted towards the red. 
Similarly, then, lines observed during (jj — 0.5 
should be blueshifted. Qualitatively this appears 
to be consistent with what is observed in Vela X- 
1 as shown in Figure 7. Quantitatively, however, 
they are inconsistent — the measured line shifts 
and widths are too small. First of all, the aver- 
age line shifts of \v\ '--^ 300 — 400 km corre- 
spond to wind velocities at approximately ^ 20% 
of the stellar radius from the photosphere, assum- 
ing a CAK profile given by Equ. 1 and ignoring 
projection effects. Second, the lines have widths 
of (7 ^ 200 km s^^ and are only marginally re- 
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Figure 9. The forest of Si K lines observed in the spectrum of Vela X-1. Ka lines from at least 8 charge 
states are detected in this narrow wavelength range. 



solved with the Chandra HEG. Finally, since most 
of the line photons observed during (jj = 0.5 come 
from a cylindrical region along the line of cen- 
ters, the magnitude of the observed velocity shifts 
should have been larger than that at (j> = 0.0, be- 
cause the wind material in this cylinder is mov- 
ing right towards the observer at = 0.5. In 
other words, the observed projected velocity of 
V ^ —200 km s~^ is probably close to the true 
wind velocity of this material. One possibility 
is that X-ray photoionization inhibits UV driv- 
ing of the wind and reduces the terminal velocity, 
in which case the wind is accelerated only up a 
certain velocity and coasts at a constant velocity 
(Hatchett & McCray 1977). Another possibihty 
is that most of the emission does, in fact, come 
from a narrow radius interval near the surface of 
the companion. In either case, the similarity in 
the line ratios between the two phases is very dif- 
ficult to understand. 

Similar puzzling behavior has also been ob- 
served in a triple star system S Ori A (09.5 II 
+ B0.5 III + early B; Miller et al. 2002) as well 
as in the isolated BO. 2 star t Sco (Cohen et al. 
2003), which show X-ray line broadening only on 



the order of cr ^ 400 km s^ . These are again a 
few times smaller than the terminal wind velocity, 
and are also smaller than those seen, for example, 
in a younger 04f star ^ Pup (Kahn et al. 2001; 
Cassinelh et al. 2001). 

Miller et al. (2002) also report evidences 
against a simple spherically-symmetric wind 
through absorption line measurements in Cyg X- 
1. Using several orbital-phase resolved spectra of 
the source obtained with Chandra, they were able 
to qualitatively infer the mass distribution in the 
system that led them to conclude that the wind 
is focused along the line of centers, most likely 
through gravitational attraction by the compact 
object. 

4. Fluorescence Line Spectroscopy 

Another very important discovery that resulted 
from the availability of high resolution data is 
the detection and identification of K-shell sili- 
con fluorescent lines from a wide range of charge 
states. Shown in Figure 9 is the HETG spectrum 
of Vela X-1 observed during (/) = 0.5. K-shell fluo- 
rescent lines from essentially all charge states are 



8 



Sato et al. 



detected in this narrow wavelength range, which 
covers only a few resolution elements of a CCD 
spectrum. This forest of lines would appear as 
one narrow (Si XIV) and one broad feature (Si II 
- XIII) as seen, for example, in the ASCA spec- 
trum of Vela X-1 (see Sako et al. 1999). 

Atomic calculations show that the strong line 
at A = 7.12 A is an unresolved blend of Ka 
lines from Si II to Si VI (fluorine- like) . The first 
few charge states (Si II and III), however, prob- 
ably do not exist in the wind of Vela X-1 due to 
the presence of a strong ionizing UV field from 
the companion star. From Si VII (oxygen-like) 
and higher, the Ka line complex shifts by ~ 70 
mA per charge state until Si XII (lithium- like) , 
and, therefore, are well-resolved with the Chan- 
dra spectrometers. 

Although detailed modeling of the Vela X-1 
spectrum is still underway (Liedahl et al. 2003), 
it is already clear that this provides an ex- 
tremely powerful spectroscopic tool for studying 
the structure of stellar winds. Since essentially 
all of the possible charge states are observed, the 
sum of the ion column densities yield the absolute 
total column density through the wind, which is 
sensitive to the total mass loss rate of the com- 
panion star. The column density distribution and 
the variability with orbital phase also allows us 
to infer the density spectrum and the clumping 
properties of the wind. 

Finally, we note that accurate elemental abun- 
dance measurements are also possible with a well- 
exposed high-resolution spectrum. As shown in 
Figure 10, for example, emission lines produced 
under similar physical conditions (neutral fluo- 
rescence, in this case) are detected from mul- 
tiple elements. Although the properties of the 
foreground as well as the embedded absorber 
may be rather complex, abundance ratios of ele- 
ments with similar atomic number can neverthe- 
less be constrained to within reasonable accuracy 
(^50%). 

5. Compton-scattered Iron Line 

The spectrum obtained during the PP observa- 
tion of GX301-2 shows a spectacular broad fea- 
ture that extends towards the redward side of the 
iron Kaline, which is identified as a Compton- 
scattered iron line in a low temperature medium. 
Compton-scattered iron lines have also been ob- 
served in the spectra of several Seyfert galaxies 
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Figure 10. Reflection-dominated spectrum of 
GX301-2 observed during the PP flare. The spec- 
trum is highly absorbed and exhibits near-neutral 
fluorescence lines from at least 7 (possibly 9) ele- 
ments. 



(references) but with much lower statistical sig- 
nificance. The spectrum of GX301-2 shown in 
Figure 11 contains ~ 6500 iron Ka counts in 
the narrow component and ~ 2000 counts in the 
Compton shoulder. 

The shape and flux of the Compton profile is 
sensitive to the electron column density and tem- 
perature of the scattering medium. The profile 
exhibits a rather sharp drop at 2Ac from the cen- 
ter of the narrow line, which implies that the tem- 
perature cannot be much higher than ~ 5 eV. 

A detailed analysis of the spectrum and its tem- 
poral variation is presented in Watanabe et al. 
(2003) using detailed Monte Carlo simulations of 
the line profile. It is found that the Ka lines 
are scattered in the same medium where they are 
produced. The similarity in the column density 
inferred from both continuum absorption and the 
Compton profile suggests that this medium sur- 
rounds the X-ray source. The covering fraction 
of this medium with respect to the neutron star 
is estimated to be fairly large (Ar2/47r ^ 1) from 
the measured equivalent widths of the emission 
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Figure 11. The iron K line observed in the spec- 
trum of GX301-2. The Ka hne clearly shows 
a broad tail towards longer wavelengths with a 
width of approximately twice the Compton wave- 
length (see inset). 



line and the depth of the corresponding K-shell 
absorption edge. 

Compton broadening can, of course, occur not 
only in emission lines, but in absorption lines as 
well. As in the emission line case, an intrinsically 
narrow absorption line produced in a Compton- 
thick medium will exhibit a broad red wing down 
to roughly twice the Compton wavelength. This 
effect might be relevant for understanding the ab- 
sorption lines observed in Cyg X-1, which show 
asymmetric wings towards the longer wavelengths 
(see, bottom panel of Figure 4). 

6. Future Work 

High resolution spectroscopic data obtained 
with the grating spectrometers onboard Chan- 
dra and XMM-Newton have provided a dynam- 
ical view of the ionized stellar wind in HMXBs. It 
is clear, however, that much work remains to be 
done. In particular, detailed physical modeling 
of each of the individual sources is likely required 
to understand the wide variety of phenomenol- 
ogy observed in the data. Radiative transfer ef- 
fects coupled with the wind dynamics may be rel- 
evant for interpreting the orbital phase variability 
of the line spectra seen in some sources. Other 
purely spectroscopic issues, such as the proper- 
ties of the Si K fluorescent lines and the Comp- 



ton recoil spectrum, must be studied in more de- 
tail to fully exploit the diagnostic capabilities of 
high resolution X-ray spectroscopy. Despite the 
different excitation mechanisms that dominate in 
HMXB systems and isolated OB stars, compar- 
ative studies of their X-ray spectra, particularly 
between similar spectral types, will be helpful for 
understanding the wind dynamics and how they 
are affected by the presence of the orbiting X- 
ray source. Although additional data will also 
certainly be helpful for understanding some of 
the questions raised with the present data, ex- 
ploratory modeling will likely provide new insight 
into the physical nature and geometry of the cir- 
cumsource media in HMXB systems. 
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